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ABSTRACT 

Active galactic nuclei with misaligned jets have been recently established as a class of high-energy gamma-ray sources. M87, a 
nearby representative of this class, shows fast TeV variability on timescales less than one day. We present calculations performed 
in the framework of the scenario in which gamma-ray flares in non-blazar active galactic nuclei are produced by a red giant or 
a gas cloud interacting with the jet. We show that both the light curve and energy spectrum of the spectacular April 2010 flare 
can be reproduced by this model, assuming that a relatively massive cloud of ^ 10^' g penetrates into the jet at few tens of 
Schwarzschild radii from the super-massive black hole. 

Subject headings: active galactic nuclei: jets - TeV photons: variability - stars: red giant 



1. INTRODUCTION 

The nearby radio galaxy M87 is a unique source for 
studies of the physics of relativistic plasma outflows and 
the conditions in the surroundings of super-massive black 
holes (SMBH). Because of its proximity (16.7 ± 0.2 Mpc; 
iMei et alJl2007l) and the very massive blac k hole at its cen- 
ter, wifli mass Mbh ^ (3 - 6) x 10^ tMacchetto et alj 
119971; iGebhardt & Thomasll2009l) . high resolution very long 
baseline interferometry (VLBI) at radio wavelengths enables 
one to directly probe structures with sizes down to < 100 
Schwarzschild radii (^sch)- From the detection of super- 
luminal features in the jet at optical and radio wavelengths, it 
has been possible to constrain the jet orientation a ngle towards 
the line of sight on the sub -kpc scale to 9 < 20° tBiretta et alj 
. [19991: lCheung"^etani2007l) . 

M87, an active galactic nuclei (AGN) with misaligned jets, 
shows very high-energy (VHE) recurrent activity with vari- 
ability timescales of a few days or less ( Aharonia n et alJ 
; l2006[ lAlbertetalJ TM. lAcciari et al.1 12009. 2010) . Re- 
I cently, a bright VHE gamma-ray flare has been simul- 
taneously detected in April 2010 by three ground based 
■ Cherenkov telescopes: H.E.S.S., MAGIC and VERITAS 
I jAbramow ski et al.ll201 II) . The detection of the VHE flare 
triggered further observations in X-rays (Chandra), and ra- 
dio (43 GHz; VLBA). The exceUent sampling of the VHE 
gamma-ray lightcurve allows a precise temporal characteriza- 
tion of the flare, which is well described by a two-sided expo- 
nential function with significantly different flux rise and decay 

times of t^^^ = 2.9 days and ■j-^^^'^^ = 0.9 days, respectively. 

The peak flux was $>().35TeV — 2 x 10~" phcm"^ s"', which 
given the unbeamed nature of the radiation allows the deriva- 
tion of a safe estimate of the luminosiy, ^ 10"^^ erg s"'. X- 
ray Chandra observations taken approximately 2-3 days after 
the peak of the VHE gamma-ray emission reveal an enhanced 
flux from the core by a factor of ^ 2 (Lx ^ 10**' erg s~'), with 
a variabihty time scale of < 2 days (lAbramowski et aLll201 lb 
lHarris etalJl2011l) . VERITAS obtained VHE spectra consis- 
tent with a power-law for three flare ph ases: rising flux, peak 
flux, and falling flux dAhu et al.ll2011h . At the peak of the 
flare, the photon index was ~ 2.2, and there is indication at 
a few a level that the spectrum is somewhat softer in the ris- 



ing and falling phases. Fermi has detected the source above 
lOOMeV with a luminosity « 5 x 10"*' erg s"' ( Abdo et aO 
l2009h . However, given the M87 flux and the Fermi sensitiv- 
ity, this instrument could not probe day-scale variability in 
this source. 

Several theoretical scenarios have been proposed to ex- 
plain the Te V flares i n M87. Generally, Georganopoulo sTt alJ 
12005) and Lenain et al. ( 2008) showed that one-zone (homo- 
geneous) leptonic synchrotron self-Compton (SSC) models 
are unlikely to explain the observed TeV spectrum of M87. 
There are also leptonic jet models with differe nt geometric 
struct ures. In the case of the multi-blob scenario jLenain et alj 
120081) a low magnetic field in the emitting region is required, 
which may be at odds with the fact th at these regions of the jet 
are Ukely strongly magnetiz ed (e.g- lKomissarov et^l.ll2007l : 
iBarkov & Komissarovl2008l) . This problem may be overcome 
if the acceleration and/or the emission processes take place in 
a weakly magnetized cloud rather than in the jefl On the 
other hand, the spine-sheath model ( Tavecchio & Ghiselliiil 
120081) predicts very strong gamma-ray absorption if it tries 
to explain the VHE hard spectra and the fast variability. 
Another possibility is the jet-in-jet model of iGiannios et alj 
(120101) , which can reproduce the spectrum of the 2010 flare, 
but does not provide at this stage a quantitati ve prediction o f 
the lightcurve of the flare. In the work by Cui et alJ (1201 ll) . 
the authors may explain the VHE flares with an external in- 
verse Compton model, but they need to assume a very wide 
jet to be able to invoke Doppler boosting. A SMBH magne- 
tospheric or igin for the TeV emission in M87 has been also 
proposed (Neronov & Aharonian '2007; Rieger & Aharonian] 
2008; Vincent & Lebohec 2010; Levinson & Rieger 2011). 
Magnetospheric models may explain a hard spectrum at VHE, 
but there is at present no detailed quantitative prediction for 
the VHE lightcu rve. Synchrotron-p roton emission may also 
operate in M87 dReimer et alj 12004 ), but to explain the ex- 
tension of the gamma-ray spectrum beyond ^ 1 TeV, strong 
Doppler boosting is needed, which is not the case in M87. 
Finally, in the paper by Barkov et al. (2010b), the TeV flares 
can be explained by the interaction between the M87 jet and 
a dense gas cloud formed out of the disrupted atmosphere of 

' Magnetic diifusion longer than the flare timescale would prevent the jet 
magnetic field from entering the cloud. 
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Accretion disc ' 

Fig. 1 . — Sketch of the considered scenario; a) penetration of an RG, with 
the external layers detached due to gravitational disruption, into the jet; b) 
penetration of a massive clump of matter into the jet. 



a red giant (RG). The emission is produced by proton-proton 
(pp) interactions between protons accelerated at the jet-cloud 
interface. This model predicts hard gamma-ray spectrum 
from GeV to TeV energies and fast variability on timescales 
of days. Unlike in models in which the VHE emission takes 
place in the dilute jet matter, pp interactions in the dense 
cloud can be energetically very efficient, and the short dy- 
namical timescale provides with the fast variability required 
to explain the flares in M87. 

In this work, we dev elop further the scenario presented 
in iBarkov et alj ( 1201 Obi) . We demonstrate that this model 
can naturally explain both the ultra short variability and the 
gamma-ray spectrum as detected. In Sect. IH we describe the 
model, the radiation features of which are explored in Sects. [3] 
and|4] and discussed in Sect.|5] 

2. THE MODEL 

The m odel considered here is based on the scenario pro- 
posed by iBarkov et alj (l2010bh . in which the envelope of an 
RG, partially tidally disrupted by the SMBH gravitational 
field, is shocked by the jet and torn away from the stellar core. 
Due to the jet impact, the RG envelope is blown up by the 
jet total pressure, forming a cloud of gas heated and acceler- 
ated downstream (see also lBosch-Ramon et al.ll20T2l) . In the 
present work, we will consider a generic gas cloud sufficiently 
massive to have a dynamical impact on the jet, at least tempo- 
rally, so it can tap a substantial fraction of the jet luminosity. 
The interaction of the jet with such a clump can convert a part 
of the jet magnetic and kinetic energy into internal energy, 
and a significant fraction of it could go to accelerate protons 
and electrons. Given the large magnet i c fields expected in the 
jet base (e.g.. lKomissarov et alj|2007l : IBarkov & Komissarovl 
I2OO8), electrons are unlikely to reach TeV emitting energies 
unless the accelerator is screened from the jet magnetic field, 
whereas protons do not suffer from this limitation. A signif- 
icant part of the accelerated protons can reach the cloud or 
the star blown material, in which case optically-thick pp col- 
lisions will lead to significant gamma-ray production in the 
early stages of the cloud expansion. A sketch of the scenario 
considered here is shown in Fig.[T] 

The atmosphere of an RG provides a good target for the 
jet to interact close to the SMBH, which is not the case for 
the stellar atmospheres of main sequence stars. The reason is 
that RGs have external layers that are much less gravitation- 
ally bounded to the stellar core. In the vicinity of a SMBH, the 



external lay ers of an RG wil l suffer significa nt tida l disruption 
(see Khokh lov et alj|1993allbl: ID iener e t 111119971: lAval etalJ 
[2000; .Ivanov et alj|2003l: iLodato et al.ii2009l) . and a mass as 
high as > 10^" g can be left almost gravitationally unbound. 
Therefore, if an RG penetrates into the innermost region of 
the jet, it can suffer the loss of its external layers due to jet ab- 
lation. Without gravitational disruption, the mass loss will 
be significantly redu ced except for very powerful jets (see 
IBarkov et alj l2010al) . Winds from stars could be effective 
clumps for the jet as well, but at larger distances from the 
SMBH, implying thus longer dynamical timescales. 

We note that, in addition to disrupted RG envelopes, other 
types of matter clumps could also be considered, like dark 
and ionized clouds near wea k and powerful acti ve galactic 
nuclei, respectively (see, e.g.. lAraudo et alj|2010l and refer- 
ences therein). One should mention in this context that a 
cloud of ^ 10^^ g wa s detected in the center of our Galaxy 
(iGillessen et al.ll2012h . 

2.1. The cloud dynamics 

One of the key parameters of the model is the power of the 
jet, which in M87 is Lj « (1 -5) x lO'*'* erg s"' (lOwen et alJ 
12000^. In this work, we fix this value to Lj « 5 x 10 erg s . 
From Lj and the jet radius, Rj = 9zj, we can derive the jet 
energy flux at the interaction height Zj'. 

Fi = « 3 X 10'^ Lj,44zr6r2 erg cm'^ s"' . (1) 

where 6'_i =6/0.1 is the jet semi-opening angle in radians, and 
2i6 = (zj/lO'^cm) is the distance from the SMBH at which the 
cloud crosses the jet. 

In the RG case, there are two tidal disruption regimes. Un- 
der strong tidal interaction, the radius of the RG, Rrg, is 
larger than t he tida l disruption radius of the star, R^,j (see 
Eq. 2 in Bark ov et al . 2010b). In that case, the RG en- 
velope becomes elongated alon g the direction of motion of 
the star (jKhokhlov et al.l[T993ah . Under weak tidal interac- 
tion, when Rrc ^ j?H,T, the envelope is still roughly spherical 
jKhokhlov et al.|[l993 b'). In both situations, the outer layers 
of the star will be swept away by the jet, forming a cloud that 
will be quickly heated up and expand. 

We study the time evolution of the cloud adopting a very 
simplified hydrodynamical model for the cloud expansion. 
The heating of the cloud is caused by the propagation of shock 
waves, which are formed by the pressure exerted by the jet 
from below. Therefore, the cloud pressure is taken similar to 
the jet pressure (regardless it is of kinetic or magnetic nature): 

F 

/^j = — ^/'c « (7-l)ec, (2) 
c 

where c is the speed of light, and 7 the adiabatic index. We 
have fixed 7 to 4/3, as corresponds to a relativistic gas, since 
the gas cloud, a t least at the initial s tages, is radiation pres- 
sure dominated (IBarkov et alj|2010b l). The cloud expands at 
its sound speed (cs), since the lateral external pressure, exerted 
by perturbed jet material, is significantly smaller than the jet 
total pr essure exerted on the cloud bottom (Pittard et al. 201()t 
iBosch- Ramon et al. 2012). As the cloud expansion proceeds, 
its pressure becomes smaller: when it is lower than the jet to- 
tal pressure, new shocks develop in the contact discontinuity, 
heating further the cloud and thus accelerating its expansion. 
We are not concerned here with the latest expansion phase at 
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which the jet lateral pressure may confine the cloud. We do 
not take into account either the energy transfer due to cosmic 
rays entering in the cloud, which would enhance the late ex- 
pansion rate. 

Numerical calculatio ns (e.g. iG regori et al.' '20001 
Nakamura et a l. 2006: P ittard etal.l20 Id: Bosch-Ramon et alj 
20121) show that the cloud is destroyed on few cloud-crossing 
times, Tco/cso, where r^o and Cso are correspondingly the 
initial radius and sound speed of the cloud after the first shock 
wave has crossed it. These simulations also show that the 
radius of the volume containing fragments of the destroyed 
cloud can grow up to an order of magnitu de compared to r^o 
(see Fig. 15 in Bosch-Ramon et aT] 120 121) . The fragmented 
cloud, with a velocity still different from that of the jet, 
continues to be suitable for shock formation and particle 
acceleration. The assumption of a spherical cloud, as adopted 
here, is a simplification, but it allows an analytical treatment 
of such a complicated system. In what follows, we present, 
depending on the tidal strength, possible analytical dynamical 
models of the cloud to be used to compute the jet-cloud 
interaction radiation. 

2.1.1. Weak tidal interaction (spherical case) 

The system of equations for the cloud radius and expansion 
velocity evolution with time, t, in the weak tidal interaction 
case, can be written as follows: 



dt 



3 3(7- l)c 
/ 7(7-l)£e ^ 

V Me 



(3) 



(4) 



where and are the cloud radius and mass, respectively. 
The solutions to Eqs. (|3]|4|i are: 
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where rco, assumed to be similar to 7?,t, is the initial cloud 
radius, and fee the cloud characteristic expansion time: 



-(. 



3cMe 
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1/2 



:5xl0^(Me28/Fy44re0,13)'^'s, (7) 



where Me28 = M^/ 10^** g. Since the velocity along the jet di- 
rection is non-relativistic during the time of interest, we ne- 
glect this motion component. 

2.1.2. Strong tidal interaction ( elongated case ) 

In the case of strong tidal interaction the RG atmosphere is 
stretched in the direction of motion of the star, and the expan- 
sion will be now cylindric; this would also apply to other types 
of clumps (see below). In such a case, the initial cloud cylin- 
drical radius r^o can be signi ficantly sma ller than the length of 
the disrupted atmosphere, (lAval et al.|[2000 ). The system of 
equations describing this case can be written as 



dt 



= Cs = 



" (7-l)c 
/ 7(7-l)£c 
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1/2 



(8) 
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Fig. 2. — Cloud radius evolution for the elongated and spherical cases com- 
bined. The parameter values are: fj = 1.6 X lO'* erg cm"^ s"', = lO" cm, 
/c = 2 X 10''' cm, and Mc = 2 X lO^' g. 



Substituting Eq. (O into Eq. (|9]), we obtain: 



(10) 



As in the weak case, r^o and fee are the initial radius and the 
expansion time of the cloud, where: 



fee — 



/ cMe 



1/2 



-1/2,-1/2 



c,28''j,14 'c,14 



'c,i4 day, 



(11) 



with Ze,i4 = (/e/10"*cm) and FjM = (fj/lO^^erg cm'^s"'). 

As mentioned above, in case the clump were not of RG ori- 
gin, the initial cloud shape is more uncertain given the lack 
of a stellar core. It seems reasonable however to assume that 
the cloud will be elongated by tidal forces and fast motion 
in the intense SMBH gravitational field. Therefore, in next 
section we adopt the dynamics of the strong tidal interac- 
tion, and when re = k we switch to the spherical solution. In 
Fig.|2l the evolution of the radius is shown for the following 
characterisic parameter values: Fj = 1.6 x 10'"* erg cm"^ s"', 
re = 10'^^ cm, Ze = 2 X lO'"* cm, and Me = 2 x 10^'' g. The value 
of Fj has been derived assuming Lj = 5 x 10'*'* erg s"', 6 = 2° 
and Zj = 3 X 10"' cm, about 20 Rsch forMeH = 6 x 10'^ Mq. 

3. RADIATION PROCESSES 

3.1. VHE gamma rays from pp collisions 

As noted in Sect.[T] it seems likely that the jet is still mag- 
netically dominated at z ^ zj . One can estimate the magnetic 
field in the jet through equipartition arguments at the level of 



8Lj 



Tj-'G, 



(12) 



where Lj 44 = Lj/lO erg s" , and Fj is the jet Lorentz factor, 
of order 1 . Such a high magnetic field, if not dissipated, will 
prevent the formation of a hydrodynamical shock . However, 
the fi -field could have an alternating polarity (e.g. iLvubarskvl 
120101) . and compression against the clump may lead to very 
effective magnetic reconnection, a potential mechanism to ac- 
celerate particles (e.g. Zenitani & Hoshino. .2001: ?). Mag- 
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netic reconnection can also lead to the eventual develop- 
ment of a hydrodynamic shock and strong turbulence (e.g. 
ISironi & Spitkovskvll201 V ). both giving rise to Fermi I and II 
type acceleration processes and to magnetic field suppression. 
Given the complexities of the processes at play, we postulate 
here that particles are accelerated in the jet-cloud interaction 
region without specifying the acceleration mechanism. 

Assuming that the magnetic field in the particle acceleration 
region is of the order of Bj, one can estimate the acceleration 
timescale: 

4cc=-^ ^~0.HE2B-\s, (13) 

where ^ is the acceleration efficiency parameter, q the parti- 
cle charge, £2 = E/10'^ TeV, and Bj.2 = Bj/IO^G, being the 
maximum energy of protons and electrons: 

Ep max « ~ 10'Bj.2rc.i4r'^' TcV (14) 

and 

£emax«,/^^«6Brf r'/'TeV, (15) 

respectively. Equation (fT4] l is obtained from limiting the pro- 
ton acceleration by Bohm diffusion escape from the interac- 
tion region, of size r^.n = (r(;/10''*cm), and Eq. ( fTSl l is ob- 
tained from limiting the electron acceleration through syn- 
chrotron cooling and neglecting other cooling factors (e.g. in- 
verse Compton). Even taking a rather high efficiency, say ^ ^ 
10 (note that for mildly relativistic shocks ^ > 100(0. 3c/v)^), 
electron energies will be too low to explain the spectrum of 
M87 measured up to energies of 10 TeV jAharonian et al.l 
l2006t lAliu et akl 1201 Ih . whereas proton energies can easily 
reach ^ 1 PeV. Even for a diffusion regime faster than Bohm 
diffusion, protons could still reach high enough energies to 
explain the obs ervations. 

As shown in iBarkov et al.l (1201 Obi) , the cloud density can 
be high, making pp interactions the best channel for gamma- 
ray production in the cloud-jet scenario. The characteristic 
cooling time for pp collisions is: 

hp- = 10^Mc.28r^HS> (16) 

CfWc 

where is the spherical cloud density and Cf is a factor that 
takes into account the cloud compression by the jet shock. 
Regarding other proton radiation mechanisms, proton syn- 
chrotron will not be effic ient due to a high-energy cutoff of 
photons around 300 GeV(Aharonian"2000), which cannot be 
Doppler boosted in the case of M87. Photo-meson production 
can be also neglected, since fp^ ^ 5 x 10^L4J ^cMCkevS, i.e. 
fp7 ^ tpp'^ where ekev = (e/1 keV) and L41 = (L/IO"*' erg s"') 
would be the energy of the target photons produced in the 
region (e.g. synchrotron X-rays; see below) and their lumi- 
nosity, respectively. 

\vi pp interactions, the fraction of the proton energy trans- 
ferred per collision to the leading 7r"-rn eson, which yields two 
photons, is « 0.17 (iKelner et al.ll2006h : this value in the opti- 
cally thick case is larger by a factor of 2. Therefore, one can 
characterize the proton-gamma-ray energy transfer by 

X = E^/Ep = 0.17[2-exp(-f App)] , (17) 
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Fig. 3. — Gamma-ray pp photon lightcurve for the weak tidal disrup- 
tion case. The parameter values are characteristic of M87: Fj = 1.6 X 
10'** erg cm"2 s"', r^o = 3 X lO'^ cm, Cf = 7, »y = 0.3, ty = lx 10** s, and 
Mc ^ 10^' g. A small pedestal flux of 2 X 10"'^ ph cm s has been 
added. 

where f can be fixed to the characteristic time of variability 
(fv), since fpp grows much faster with time. 

Two phases of the cloud expansion can be distinguished: 
the radiatively efficient regime, i.e. with x ~ 0.34 or tpp < fv, 
and the radiatively inefficient regime, with x = 0. 17 or tpp > ty. 
Thus, from the simplifications above, the gamma-ray lumi- 
nosity in the pp optically-thick case can be written as: 

w0.34r;7rr^Vj, (18) 

where it has been assumed that the injected proton luminosity, 
Lp, grows like the jet-cloud interaction surface, oc r^; 77 is the 
efficiency of energy conversion from the jet to the accelerated 
particles. In the optically-thin regime, only a fraction ^ t^/tpp 
of Lp is lost through pp collisions, and cx . The gamma- 
ray luminosity during the cloud-jet interaction is: 

L^KiTTr]xrlFj(l- e-'^l'"' ) . (19) 

Given the fast expansion of the cloud, either in the spheri- 
cal or the elongated case, one can expect a sharp spike in the 
lightcurve. 

Analytical lig htcurves similar to that derived in 
iBarkov et alj (12010b) are presented in Figs. [3] and |4] for the 
weak and the strong tidal interaction case, respectively. The 
adopted parameter values in Fig. [3]are: Cf = 7, fv = 7 x lO** s, 
Mc PS 10^"^ g, and = 1.6 x lO''* erg cm"^ s"'. As seen from 

the figure, the VHE peak is reached at fpeak ~ 4 x 10^ s, 
with a width of ^ 1 - 2 days, adopting 77 = 0.3. To obtain a 
short rise phase we use a relatively large cloud initial radius: 
= 3 X lO'-' cm. In the strong tidal interaction case one 
can set a smaller initial radius of the RG, since the clump is 
already elongated due to tidal forces before penetrating into 
the jet. The lightcurve for the strong tidal interaction case is 
shown in Fig. 2) The adopted parameter values are the same 
as in Fig. [3]but fv, fixed now to 9 x 10"* s, and = lO'-' cm 
and Zc = 2 X lO''* cm. As seen in the figure, the VHE peak is 
reached at fpeak ~ 5 x 10^ s, with a width of ^ 1-2 days. 

3.2. Non-thermal electrons and positrons 
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Fig. 4. — The same as in Fig.|3]but for the strong (elongated) tidal disrup- 
tion case and taking ?v = 9 X lO^s, r^o = lO" cm, and Ic = 2x lO'* cm. 

Secondary electrons and positrons (e^) are injected into the 
cloud through the decay of charged secondary vr-mesons pro- 
duced in pp collisi ons. The energ y rate of the injected 
pan's IS w /2 (Kelner et al.ll2006l) . The region is quite com- 
pact, and depending on the cloud magnetic field most of the 
secondary emission can be synchrotron or SSC, with a mi- 
nor contribution from relativistic Bremsstrahlung. The initial 
magnetic field in the cloud can be relatively small; e.g. in the 
case of a (non-disrupted) RG atmosphere, the magnetic field 
is of several gauss (e.g. Ko nstantinova-Antova et al] 120081 
|20TO|) . The quick expansion can decrease the Z?-strength 
rapidly, although the continuous pumping of jet energy, plus 
complex dynamo effects in the cloud, may prevent the B- 
field from decreasing, and may even enhance it. Assum- 
ing that X-rays right after the peak come from the cloud, 
an estimate of the cloud magnetic field can be done assum- 
ing that synchrotron losses dominate over SSC. This yields 
a lower limit for the cloud magnetic field (Be), since its en- 
ergy density, B^/Stt, is to be larger than the synchrotron one, 
CDsyn S i'7X /47r r^c, where L^x is the gamma-ray luminosity at 
the X-ray observation (see, e.g., Fig.O. This yields a (loose) 
constraint on B^: 

B, ~ (8^C3sy„) - 30Lill,^r:\^G, (20) 

where L^xai = (L^x/^ x 10"*' erg s"'). Note that even for 
lower Bc-values, gamma-rays could be absorbed through 
photon-photon (77) collisions in the synchrotron field. Given 
the different possible Bc-values and radiation outcomes, we 
will treat Be after the cloud expansion as a free parameter con- 
fined to a value range ^ 0.01 - 100 G. 

Primary pairs carried by the jet may be also acceler- 
ated together with protons at the jet-cloud contact discontinu- 
ity. Since their losses are expected to be synchrotron domi- 
nated given the high Bj-value, their emission should not over- 
come the X-ray fluxes detected few days after the VHE peak. 
The primary synchrotron radiation should not absorb the VHE 
photons either All this sets an upper limit in the primary 
injection rate at ~ 0.01 around the VHE maximum. Be- 
side this restriction, however, the primary leptonic population 
is rather unconstrained, whereas secondary pairs are almost 
completely constrained by the relativistic proton spectrum and 
estimated cloud conditions, with only the Bc-value remain- 



ing free. Below we examine under which conditions these 
secondaries do not lead to strong gamma-ray absorption, and 
explore whether they may still explain the enhanced X-ray 
emission detected few days after the VHE maximum. 

4. MODELING THE HIGH-ENERGY EMISSION 

The emitting proton population has been modelled adopt- 
ing a spatially homogeneous (one-zone) emitter, in which rel- 
ativistic protons are injected with power-law energy distribu- 
tions of two types: Qp{E) oc E'"^ plus an exponential cutoff 
at /ic = 200 TeV and a low-energy cutoff at Eue = 1 TeV; and 
Qp{E) oc E~^'^ plus an exponential cutoff at E^ = 50 TeV. The 
first case would correspond to a fairly standard phenomeno- 
logical assumption plus a low-energy cutoff, and the second 
one may be associated to non-linear shock acceleration with a 
large compression ratio. 

The injection luminosity of protons has been taken Lp = 
r/TT rdt)^ Lj, and thus the model is inhomogeneous in time. 
To derive the proton energy distribution for different f-values 
(N(E,t)), the time-dependent proton injection and cloud con- 
ditions have been modelled as follows. Protons injected at a 
certain time f,, Q'p(E), are evolved for a (5f,- under the corre- 
sponding cloud conditions. The evolved population is added 
to N^^.^y(E), i.e. the accumulated population from t[ to f,_i 

(^p7ev(^)) evolved also for a St,. This numerical technique 
provides a correct solution at any time of interest, provided 
that Sti <C where t'^^^^ is the dominant evolution timescale 
of protons at f,-. At the relevant time interval, the proton evolu- 
tion is dominated hy pp interactions over other radiation cool- 
ing mechanisms. Adiabatic losses, given the roughly constant 
cloud pressure during the relevant expansion phase, do not 
play a significant role. 

In Fig. |5] we show the VHE lightcurve computed using 
N(E,t), derived adopting the Qp{E) oc E~^ case, and a proton- 
to-gamma-ray energy fraction per collision of 0.17. The re- 
maining parameter values are the same as in Fig. |2]plus a 
non-thermal efficiency r/ = 0.5. As seen in the figure, the com- 
puted VHE peak can be slightly broader than the observed 
one, although a stronger decay can be naturally explained by 
the cloud escaping the jet and shutting off the proton injection. 
The Keplerian speed at 3 x 10'^ cm (Mbh = (3-6) x IO'^Mq) 
is w 4-5 X 10^ cm s"\ and a cloud with the adopted mass 
will leave the jet after several days for Rj = lO' ^ cm. Lighter 
clouds would be dragged downs tream the jet (i Barkov et al.l 
1201 Oat iBosch-Ramon et ani2012b . To model the effect of the 
cloud escape, two possibilities have been considered together 
with cutting off the proton injection: either the cloud expan- 
sion continues at the same (growing) rate, or it is kept con- 
stant after the escape. The VHE spectral energy distributions 
(SED) obtained adopting Qp{E) oc E~^ and E~^-^ are presented 
in Fig. |6] Despite the hard proton spectrum, the lower energy 
exponential cutoff of the case with Qp(E) oc E~^'^ renders a 
very similar pp SED at VHE. 

Some remarks should be done regarding the energetics and 
the physical consistency of the discussed dynamical-radiation 
model. For either the Qp{E) oc E~^ or the Qp{E) oc /i"'^ case, 
the required Lp -value to explain the VHE event in the con- 
text of the adopted dynamical model amount^ ^ 10'*^ erg s~' . 
This means that the average efficiency of proton-to-photon en- 

^ For the case with Qp(E) oc E^^ with a low-energy cutoff at ~ GeV ener- 
gies, the luminosity budget would roughly double. 
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M^=2xl(/' g, F -1.6xI0'* erg/em's. r|NT="-5' R cm. 1-2x10" cm 




1 IdaysJ 



Fig. 5. — The VHE-luminosity lightcurve of the flare assuming proton in- 
jection all the way up to 12 days (thin dashed line) and only during 7.3 days 
(dotted -increasing radius growth speed- and thick solid line -constant radius 
growth speed-), as expected if the cloud leaves the jet. The parameter values 
are the same as in F ig.[2]plus ri = 0.5 and Qp( E) oc E^^. The data points have 
been obtained from lAbramowski et alj j2011l) . 



M.=2xllf ' g. Lpl.fixio" erg/cm\, 11^.^=0.5. R^=lo'"' cm, l^=2xl()'* cm 




log ( E ITeV] ) 



Fig . 6 . — The computed SED for pp gamma rays at the VHE maximum for 
two different injection spectra: Qp(E) oc E"^ (£le = 1 TeV, iic = 200 TeV; 
solid line); and Qp(E) oc E^''^ (no Ele, Ec = 50 TeV; dot-dashed line). The 
remaining param eter values are the same as in Fig. |3] The data points are 
obtained from lAliu et alj <2011D . 



ergy transfer is ^ 15%. This also implies a rather tight energy 
budget, ij = 0.5, although this is partially due to the limitations 
of the analytical dynamical model. For instance, the compres- 
sion ratio has been fixed to Cf = 7. However, in reality, further 
expansion makes the cloud optically thin. Since the cloud 
shocks are radiative, the compression ratio then can be po- 
tentially much higher than 7. Another effect to consider here 
is the cloud fragmentation, which can take place at the late r 
stages of the cloud expansion (see Bosc h-Ramon et aO2012l) . 
Cloud fragmentation effectively implies a larger cloud cross 

section by a factor ^ Nf^^, where A'f is the number of cloud 
fragments. Because of these facts, the real energetic require- 
ments may be relaxed by a factor of several. Note that, for 



instance, the very high and sharpe peak observed by H.E.S.S. 
could be explained when accounting for the aforementioned 
stronger compression and fragmentation. We did not include 
these effects in the present calculations since this would imply 
additional free parameters. It is also noteworthy that the back 
reaction of accelerated particles, which can heat the cloud, 
accelerate its expansion, and shorten the VHE emission rise 
time, has been neglected. To take into account all these ef- 
fects, a more accurate treatment is needed that would provide 
the coupling between acceleration, radiation and magnetohy- 
drodynamics (MHD). This will be left for future studies. 

4.1. X-rays 

Few days after the maximum in the VHE lightcurve, Chan- 
dra observed the M87 core and found it in a high-flux state, 
which decreased by a factor of two in the following days 
(|Harris et al. 2011). Although the X-ray emission does not 
necessarily come from pp-e^ pairs (plus a 77 -e^ contri- 
bution; see below), and may have a primary leptonic origin, 
it is interesting to explore whether pp interactions alone can 
explain the X-ray and VHE observations. Therefore, together 
with the pp gamma rays, we have also calculated the broad- 
band emission of the pp-e^ pairs. For that, we derived firs t 
the injection spectrum Qg±(E) using .Kelner et all (l2006l) . 
and computed its evolution at the cloud conditions around the 
VHE maximum. 

The dominant cooling channels within the cloud are syn- 
chrotron or SSC, depending on B^. Other cooling channels, 
like relativistic Bremsstrahlung or inverse Compton with IR 
photons from the M87 nucleus, are less important. Since the 
leptonic cooling processes are quite fast at the relevant en- 
ergies, the cloud radius and pair injection can be consid- 
ered as roughly constant. Under SSC dominance, however, 
the emitter still requires an inhomogeneous treatment in time 
because the cooling of pairs depends on their own syn- 
chrotron radiation, making the evolution non-linear. To ac- 
count for this, we have applied to the pairs the same nu- 
merical technique used for protons. Note that to compute the 
whole evolution of the secondary pairs and their emission, 
the injection should change in time as the pp gamma-ray 
luminosity does. As it is computationally rather expensive, 
we have restricted ourselves to the leptonic SED for the cloud 
conditions at the VHE maximum, and adopted two Bc-values 
to obtain SEDs for the synchrotron and the SSC dominated 
cases. 

As discussed in Sect. [3] effective secondary synchrotron 
emission, high enough to explain Chandra April 11th 2010 
data, can suppress the VHE emission. This can happen even 
for a very hard injected population, since synchrotron 
cooling will lead to iv oc z/""^ below X-rays. In addition, 
since a relatively high Be -value is required for synchrotron ra- 
diation to overcome SSC, electromagnetic cascades will not 
be able to alleviate the strong gamma-ray absorption. The 
situation is actually complex, since in fact the 77 pairs can 
dominate over the pairs from pp interactions, yielding even 
brighter and harder X-rays. We find that Z?c < 01 G is neces- 
sary to avoid strong absorption, and Be ^ 10 G could explain 
the X-ray fluxes few days after the VHE maximum. 

The computed SED for a low Bc-value (0.03 G), able to 
reproduce the observed VHE SED avoiding significant ab- 
sorption, is presented in Fig. [T] (top) for protons following a 
Qp{E) oc E~^. The X-ray high state can be roughly explained 
by an increase in B^, which may happen via mixing of the 
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strongly magnetized jet material with the cloud when it frag- 
ments, or by dynamo effects in the complex MHD flow. If 
this were the case, the decay of the VHE lightcurve could 
not be explained by jet escape, since high X-ray fluxes are 
only possible under sustained high pp collision rates. This 
decay though can be naturally explained by 77 absorption. 
The high Be case is shown in Fig. [T] (middle), for which we 
have adopted = 30 G. The 77 - synchrotron contribu- 
tion is also included in the figure. The SSC and relativistic 
Bremsstrahlung levels, not shown, are similar to those of the 
pp-e^ pairs. The multiwavelength SED for protons with 
Qp{E) oc E"'^ and Be = 0.03 G is presented in Fig. |7] (bot- 
tom), and is similar to that of the Qp(E) oc case (the same 
applies for the case with = 30 G, not shown). 

To illustrate further the impact of Be on gamma-ray sup- 
pression and the cloud (secondary) leptonic population, we 
present in Figs. |8] and |9] the opacities to gamma-rays, and the 
SED {E^N{E)) of the injected pp- and 77-e^ pairs, for the 
two Be-values adopted. It is worth noting that 77 absorption 
softens the VHE spectrum in Fig. |7] (middl e), in g ood agree- 
ment with observations in the decay epoch iAliu et al. 201 1^ 
Note also that the non-absorbed VHE lightcurve peak is broad 
enough to accommodate the two days separation between the 
VHE maximum and the X-ray observation. 

4.2. Thermal emission from the cloud 

As shown in lBarkov et al.l(l2010bh . the shocked cloud is ini- 
tially optically thick to its own radiation. Given the cloud 
pressure ^ very high UV photon densities will completely 
quench all the VHE radiation. With expansion, the cloud 
radiation field quickly dilutes and heats up, with a fast de- 
crease of the 77 opacity. At the peak of the flare, the cloud 
is fully ionized, moderately optically thin (je-y ~ 0.6 re_i4ne.io; 
ic,io = ne/10'"cm~^), and emitting through free-free radiation 
with temperature ^ lO"' K and luminosities of ^ 10^' erg s"'. 
The optical depth fo r TeV photons can be estimated as (see 
iBarkov etaO2010bl for further details): 



a) 



^7 7, therm ' 



10 ■'''c, 14 "c, 10^10 



1/2 



(21) 



The opacity due to thermal radiation can be therefore ne- 
glected, and the thermal radiation itself becomes important at 
MeV energies around the VHE maximum, when the jet-cloud 
energy transfer is most efficient. 

5. DISCUSSION 

The interaction of a gas cloud, or the atmosphere of a dis- 
rupted RG, with the base of an AGN jet leads to the forma- 
tion of an interaction region. There, jet energy can be dissi- 
pated in the form of relativistic protons, which can penetrate 
into the cloud, initially dense enough to render pp interac- 
tions efficient. This yields gamma rays and other secondary 
particles, in particular pairs. The energy transfer into 
the cloud, through MHD interactions (and injection of cos- 
mic rays), leads to a quick expansion that increases the cov- 
ered section of the jet and thereby the gamma-ray emission. 
Around the point when pp collisions become optically-thin, 
the gamma rays reach their maximum after a quite sharp rise. 
After that moment, pp collisions become strongly inefficient 
quenching the gamma-ray emission. The drop of the gamma- 
ray flux can be even more abrupt because of the cloud leav- 
ing the jet, and/or 77 absorption. Adopting Me ~ 10^' g and 
Lj w 5 X 10'*'* erg s~', the model can reproduce rather nicely 
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Fig. 7. — a): The computed SED for secondai^ synchrotron and SSC. 
and pp gamma rays, adopting Qp(E) oc E^^, £le = 1 TeV, Ec = 200 TeV. and 
Be = 0.03 G. The other parameters are the same as those adopted to obtain 
Fig.|6] b): The same as in a) but for = 30 G, plus the -yy—e^ synchi'otron 
emission. The approximate level of observed X-ray emission in April 11th 
2010 is also shown. The computed synchi'otron emission from the 77 -e^ 
pairs are also shown, c): The same as in a) but for a proton injection spectrum 
oc E-' ^, no Ele, and Ec = 50 TeV. 
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Fig. 8. — The computed 77 opacities for the same parameter values of 
Fig.[6] and Be = 0.03 (solid hne) and 30 G (dashed line). 
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Fig. 9. — The computed SED of the injected pp— and 77 -e^ pairs for the 
same parameter values of Fig. |6](Qp(£) oc and Be = 0.03 (solid line) 
and 30 G (dashed line). 



the lightcurve and the spectrum of the VHE flare detected in 
M87 in April 2010. 

At the VHE peak, a low Be is required to avoid gamma- 
ray absorption in the secondary synchrotron field. This is not 



problematic, since the cloud may be initially weakly magne- 
tized. The enhanced X-ray flux few days after the VHE maxi- 
mum is consistent with our model if Be increases, e.g. through 
cloud-jet mixing or complex MHD processes. A strong pre- 
diction of the proposed model is that, at the highest VHE 
fluxes of an event like the April 2010 one, to avoid absorp- 
tion X-rays should not be significantly enhanced regardless 
their origin, either primary or secondary leptons. 

The model presented here cannot be applied for particles 
whose energy evolution timescale is either longer, or strongly 
coupled to the dynamical timescale of the cloud, which pre- 
vents to derive predictions in radio. However, the compact- 
ness of the source already shows that the radio emission will 
be likely self-absorbed during the flare. In the IR/optical, 
the computed fluxes are below the observed ones even for a 
high Bc-field. Later, primary and secondary leptons may non- 
negligibly contribute to the radio and IR/optical emission, but 
a proper account of these radiation components would require 
a very detailed model of the dynamical and radiative proper- 
ties of the interacting flows. In any case, in what concerns 
the April 2010 flare, the radio and optica l emission did not 
seem to correlate with the VHE emission (lAbramowski et all 

unni). 

The hard proton spectrum adopted implies rather modest 
GeV fluxes, although if the injection proton spectrum cx 
went down to ^ 1 GeV in proton energy, the 0.1-10 GeV lu- 
minosity would be a factor of a few times higher than the year- 
averaged flux found by Fermi. However, for these fluxes and 
the sensitivity of this instrument, it seems difficult to probe 
day-scale variability. 

It is important to emphasize that, because of the lack of 
beaming and the limited jet luminosities in local Universe 
AGN, single interactions of clumps with jets i n misaligned jet 
sources (for blazars, see Barkov et al.ll2010al) are hard to de- 
tect. Only very nearby objects, like M87, Cen A, and probably 
NGC1275, could provide detectable radiation of this origin. 
More distant though still local objects may be also detectable, 
since more clumps are available at further distances from the 
SMBH. This emission would appear as persistent though (e.g. 
lAraudo eTa l. 2010; Bosch-Ramon et al. 2012), and the mod- 
est jet power of the potential candidates would probably re- 
quire rather long observations. 
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